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ABSTRACT 

We study the coohng of superfluid neutron stars whose cores consist of nucleon matter 
with the Akmal-Pandharipande-Ravenhall equation of state. This equation of state 
opens the powerful direct Urea process of neutrino emission in the interior of most 
massive neutron stars. Extending our previous studies (Gusakov et al. 2004a, Kaminker 
et al. 2005), we employ phenomenological density-dependent critical temperatures 
Tcp(p) of strong singlet-state proton pairing (with the maximum T™'^^ ~ 7 x 10^ K 
in the outer stellar core) and Tcnt(p) of moderate triplet-state neutron pairing (with 
the maximum ~ 6 x 10^ K in the inner core). Choosing properly the position of 

T^f^ we can obtain a representative class of massive neutron stars whose cooling is 
intermediate between the cooling enhanced by the neutrino emission due to Cooper 
pairing of neutrons in the absence of the direct Urea process and the very fast cooling 
provided by the direct Urea process non-suppressed by superfluidity. 

Key words: Stars: neutron - dense matter 



1 INTRODUCTION 

Fundamental properties of supranuclear matter in the cores 
of neutron stars, such as the equation of state and critical 
temperatures of baryon superfluidity, a re still poorly known. 
In pa rticular, they can be studied (e.g.. lYakovlev &: Pethic^ 
I2OO4I and references therein) by comparing the theory of 
cooling isolated neutron stars with observations of therma l 
radiation (e.g.. IPavlov et aljliool IPavlov fc ZavlinI l2003l) 
from these stars. 

In this paper we study the cooling of superfluid neu- 
tron stars with the cores composed of neutrons, protons, 
electrons, and muons. We employ model equations of state 
(EOSs) of neutron star cores based on the EOS of Akmal, 
Pandharipande & Ravenhall (1998) (APR, Section 4) which 
is currently assumed to be the most elaborated EOS. It 
opens the powerful di rect Urea process of neutrino emission 
iLattimer et alj[l99 j) in neutron stars with masses slightly 
lower than the maximum mass. Our cooling models will be 
based on standard (nucleon) composition of neutron stars 
and standard neutrino processes (no "exotic" physics in- 
volved) . 

T he p aper extends our previ ous studies iGusakov et alJ 
l2004al and lKaminker et al.ll2005l : hereafter Paper I and Pa- 
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per II, respectively) of cooling neutron stars with superfluid 
nucleon cores, where the direct Urea process is forbidden. In 
this scenario, observations of cold neutron stars (e.g., of the 
Vela pulsar) are explained by the enhanced neutrino emis- 
sion due to Coope r pairing of neutrons in neutron star cores. 
iPage et al.l (jio^l), who suggested a similar cooling scenario, 
called it "Minimal Cooling". Note, however, that the real- 
ization of this scheme by Page et al. and in Papers I and II 
is different (see Paper II for details). 

In Papers I and II we employed the EOS of 
iDouchin fc HaenseJ J200j) (the DH EOS, which forbids the 
direct Urea process in all stars) and a model of triplet-state 
neutron pairing with the maximum of the critical tempera- 
ture 6 X 10* K in the inner stellar core. This pair- 
ing is needed to produce an enhanced neutrino cooling of 
massive stars and explain the observations of neutron stars 
coldest for their age (such as the Vela pulsar). In Paper I 
we assumed also strong singlet-state proton pairing in the 
outer stellar core, with T™^" > 5 X 10*^ K. It slows down 
the cooling of low-mass neutron stars and enables us to in- 
terpret the observations of stars warmest for their age (e.g., 
RX J0822-4300 or PSR B 1055-52). In Paper II we included 
the effects of surface (accreted) envelopes of light elements. 
They allowed us to explain the observations of warmest stars 
assuming weaker proton pairing {T^^^ > 10^ K). 

Now we will add the effect of the direct Urea process in 
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most massive stars in accordance with the APR EOS. We 
will use similar (but slightly different) models for strong pro- 
ton pairing and moderate neutron pairing in the stellar core. 
These superfluid models are sufficient to explain current ob- 
servations of all cooling neutron stars without opening the 
direct Urea process. However, most massive stars predicted 
by the APR EOS can be much colder than the stars observed 
today and, we will focus on their properties. For simplicity, 
we will neglect the effects of accreted envelopes. They are 
unimportant for massive neutron stars; we could easily in- 
clude them as in Paper II. 



2 OBSERVATIONS 

Table f summarizes the observations of isolated (cooling) 
middle-aged neutron stars (10^ yr < t < 10'' yr), whose 
thermal surface radiation has been detected or constrained. 
We present the estimated stellar ages t and the effective 
surface temperatures (redshifted for a distant observer). 
The data are described in Paper I and slightly corrected in 
Paper II. 

The surface temperatures of some sources (labeled by 
"') have been obtained by fitting their thermal radiation 
spectra with hydrogen atmosphere models. Such models are 
more consistent with other information on these sources 
fe.g.. lPavlov fc Zavlin 2003.) than the blackbody model. For 
other sources (e.g., for Geminga and PSR B1055-52, labeled 
by ) we present the values of inferred using the black- 
body spectrum because this spectrum is more consistent for 
these sources. The surface temperature of RX J1856. 4-3754 
is still uncertain. Following Paper I we adopt the upper limit 

< 0.65 MK. Finally , T,°° for RX JQ720.4-3125 is taken 
from lMotch et all i2003r) . who interpreted the observed spec- 
trum with a model of a hydroge n atmosphere of f inite depth. 

Note also new results by iKargaltsev et all ll2005il for 
Geminga presented in Table 1. As in Papers I and II we 
retain 20% errorbars for this source. 



3 COOLING WITHOUT THE DIRECT URCA 
PROCESS 

We will simulate neut ron star co oling using our general rel- 
ativistic cooling code iGnedin et al. 2001^') which calculates 
cooling curves (the dependence of on the stellar age t). 
At the initial cooling stage {t < 100 yr) a star cools via neu- 
trino emission but the stellar interior is non-isothermal. At 
the next stage (10^ yr < t < 10^ yr) the neutrino emission 
dominates but the interior is isothermal. Later {t > 10^ yr) 
the star cools mainly via the surface photon emission. 

For the adopted EOSs, all constituents of dense mat- 
ter (nucleons, electrons, and muons) exist everywhere in the 
stellar core, except for muons which appear at p ~ 2.0 x 10^* 
g cm~'^. As in Paper I we use the relation between the 
effective surface temperature and the temperature at the 
bottom of the heat-blanketing iron envelope calculated by 
[Potekhin et aL (.200^1 . 

Neutron star cooling is strongly affected by superflu- 
idity of nucleons in stellar interiors. Any superfluidity is 
characterized by a density-dependent profile of critical tem- 
perature Tc{p). In principle, an EOS and a superfluidity 



model, employed in cooling simulations, have to be obtained 
from the same nuclear Hamiltonian. However, such calcula- 
tions are currently absent and one usually treats the EOS 
and superfluid properties as independent. Although this ap- 
proach is simplified it sounds reasonable at the present 
stage of the investigation because the EOS is the bulk 
property of the matter (determined by the entire Fermi 
see of baryons) while superfluidity is a Fermi-surface phe- 
nomenon and can be, therefore, relatively independent. Mi- 
croscopic theories predict two main types of superfluidity 
inside the neutron star core: singlet-state proton pairing 
(Tc = Tcp) and triplet-state neutron pairing {Tc = Tent). 
These theories give a large scatter of critical tempera- 
tures depending on a nucleon-nucleon interaction model and 
a ma n y-body theory e mploy ed (e.g., iLombardo fc Schulzd 
2001', 'Yakovl ev efaU \l99^ see also recent papers by 
Schwcnk & F rirnanI |2004 iTa katsuka & Tamagaki 20o|, 
Zuo et all2004lTanigawa et alhooi) . As in Papers I and II, 
we will treat T'cp(p) and rcnt(p) as phenomenological func- 
tions of p which do not contradict some microscopic calcu- 
lations. These functions could be (in principle) constrained 
by comparing theoretical cooling curves with observations. 

Superfluidity of nucleons affects the nucleon heat capac- 
ity and suppresses neutrino processes involving superfluid 
nucleons: Urea and nucle on-nucleon bremsstra hlung pro- 
cesses (as reviewed, e.g., bv lYakovlev et al Moreover, 
superfluidity induces an additional neutrino emission mech- 
anism associated with Cooper pairing of nucleons (Flowers, 
Ruderman & Sutherland 1976) . All these effects of superflu- 
idity are included into our cooling code in the same manner 
as in Papers I and II. 

In Fig. we plot two models of nucleon pairing em- 
ployed in the present paper, model p for strong singlet-state 
proton pairing (with the peak of Tcp (p) approximately equal 
to r™"'' f» 6.8 X 10^ K); and model nt for moderate triplet- 
state neutron pairing (T^t" ~ 6 x 10* K). Models p and 
nt slightly differ from similar models pi and ntl in Pa- 
pers I and II. Note that the curves Tcp(p) and rcnt(p) are 
almost insensitive to the employed EOS: DH (Papers I and 
II) and APR I, II, and HI (Section 4). Our pairing model 
nt seems specific (shifted to too high densities p at which 
neutron pairing dies out according to many microscopic cal- 
culations). However similar models have been obtained from 
some microsc opic theorie s (e.g., sec the curve m* — 0.73 in 
Fig. 1 of .Takatsuka fc Tamagaki 1997') . 

Figure 2 illustrates the cooling scenario of Paper I 
(to be compared with new scenarios). We plot a family 
of cooling curves calculated for the DH EOS, which for- 
bids the direct Urea process in all stable neutron stars 
(M ^ Mmax = 2.O5M0). The curves refer to neutron stars 
with (gravitational) masses M from 1.00 to 2.04 Mq (from 
top to bottom). 

The curves explain all the data from Table 1. The cold- 
est objects can be interpreted as most massive neutron stars 
with moderate neutron pairing in the inner cores (at densi- 
ties p ^ 8 X 10^* g cm~^). When the internal temperature 
of a massive star falls below T^^^, the cooling is noticeably 
enhanced by the neutrino emission due to Cooper pairing 
of neutrons. This enhancement is sufficient to explain the 
observations of neutron stars coldest for their age, partic- 
ularly, PSR J0205-h6449, RX J0007.0-h7303, the Vela and 
Geminga pulsars. Note that the data on RX J0007.0-(-7303 
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Figure 1. Density dependence of the critical temperatures for model p of singlet-state proton pairing and model nt of triplet-state 
neutron pairing in a neutron star core. Three vertical dot-and-short-dashed lines indicate the thresholds of the direct Urea process for 
the EOSs APR I, II, and III discussed in Section 4. The vertical dot-and-long-dashed line shows the central density of a maximum-mass 
neutron star for APR I; two other maximum-mass lines are almost the same. 
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Figure 2. Observational limits on surface temperatures of twelve isolated neutron stars (Table 1) compared with theoretical cooling 
curves calculated for the DH EOS (no direct Urea process) and for models p and nt of nucleon pairing (Fig.QJ. Fifty three curves refer 
to neutron stars of masses from 1.00 to 2.04 Mq (from top to bottom, with the step 0.02 Mq). 
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Table 1. Observational limits on surface temperatures of isolated neutron stars 



Source 



t [kyr 



[MK] 



Confid. of T„°° References 



PSR J0205+6449 (in 3C 58) 0.82 <1.02 ^) 99.8% 

Crab 1 <2.0 99.8% 

RX J0822-4300 2-5 1.6-1.9 90% 

IE 1207.4-5209 3-20 1.4-1.9 90% 
RX J0007.0+7303 (in CTA 1) 10-30 < 0.66 

Vela 11-25 0.65-0.71 68% 

PSR B1706-44 ~17 0.82lg;°^ 68% 

PSR J0538+2817 30 ± 4 ~ 0.87 

Geminga ~340 ~ 0.5 

RX J1856.4-3754 ~500 <0.65 

PSR B1055-52 ~540 ~0.75 

RX J0720.4-3125 - 1300 ~ 0.51 



Slane et al. C2004') 
Weisskopf et al. (2004) 
Zavlin et al. (1999") 
Za vlin et al. (2004) 
Halpern et al. (2004) 
Pavlov et al. C20021 
McGowan et al. (2004) 
Zavlin fc Pavlov (2004) 



Kargaltsev et al. (20051 
see Paper I 

Pavlov fc Zavlin (2003) 
Motch et al. (2003) 



^) Inferred using a hydrogen atmosphere model 
Inferred using the black-body spectrum 



are explained in this way only marginally. However, tak- 
ing into account that real uncertainties of Tg°° and t can 
be larger than the adopted ones we can accept this ex- 
planation as successful. All stars with the central density 
higher than the density at which neutron pairing dies out 
(Pc ^ 2 X 10^^ g cm~'^ in Fig. cool nearly as fast as the 
star with M = Mmsoc- This gives the lower dense bunch of 
cooling curves of massive neutron stars (with pc >2 x 10^^ 
g cm" 3) in Fig.H 

There is the other (upper) dense bunch of cooling curves 
of low-mass stars (M < 1.1 M0). The central densities of 
such stars are pc ^ 8 x 10^* g cm""^. Proton pairing is strong 
(Tcp(/5) > 3 X 10^ K) everywhere in their cores (FigEJ. It 
completely suppresses all neutrino processes involving pro- 
tons and substantially slows down the cooling. 

Finally, there is a class of medium-mass neutron stars 
which show intermediate cooling. Their cooling curves fill 
in the space between the upper dense bunch of cooling 
curves (low-mass stars) and the lower dense bunch (high- 
mass stars). We can treat PSR B1706-44, PSR J0538+2817, 
and RX J1856. 4-3754 as medium-mass stars. 

Although this scenario is consistent with all current ob- 
servations, colder stars could be detected in the future. An 
indirect evidence for the existence of such stars is given by a 
non-detect ion of compact cen tral objects in four supernova 
remnants an et al.ll2004^ . If neutron stars are present 

there, they should be really cold. It would be impossible to 
explain them with the cooling scenario presented in Fig. 2, 
but they could be interpreted adopting the APR EOS. 



4 COOLING OF APR NEUTRON STARS 
4.1 Equations of state 

Let us employ three EOSs of dense matter, APR I, APR II, 
and APR III. They represent a useful parameterization 
of the well-known EOS obtained by IXkm al et al. (1991) 
(model Argon VI S + 5v + UIX*). This parameteriz ation 
was proposed by iHeiselberg fc Hiorth-Jensenl (ll999D . We 
use the same parameters for the compressional modulus 
of the nucleon matter and vary only the parameter 7 in 



the expression for symmetry energy, S{nh) = So (rib/no)'', 
where So = 32 MeV; rib is the baryon number density, and 
no — 0.16 fm~"^ is the baryon number density in saturated 
nuclear matter. Specifically, we have taken 7=0.6, 0.575, and 
0.643 for the APR EOSs I, II, and III, respectively. These 
EOSs differ mainly in the threshold densities for opening the 
direct Urea process (see Fig. 1). 

Some properties of the EOSs and associated neutron 
star models are illustrated in Figs. and |3 and in Table 
2. For comparison, in Fig. El and Table 2 we include also 
the data on the DH E OS IPouchiii fc Haenseil200ll) and 
the original APR EOS jAkmal et al.lll998^ . The APR EOS 
corresponds to the nucleon matter with the appearance of 
neutral pion condensation. To plot th is EOS in Fig. El we 
use th e parameterization proposed by iHaensel fc PotekhinI 
(|20^), with the values of fit parameters provided by A.Y. 
Potekhin (2004, private communication). 

Table 2 gives masses M, central densities pc, and cir- 
cumferential radii 71 of two stellar configurations. The first 
configuration is the most massive stable neutron star. The 
values of A/max indicate that all five EOSs are moderately 
stifle and similar. The second configuration refers to the onset 
of the direct Urea process, (pc — po, M = Mu); it is absent 
for the DH EOS. The APR III EOS imphes larger symmetry 
energy at supranuclear densities (larger 7) and larger pro- 
ton fraction than the EOSs APR I and II. Accordingly, the 
direct Urea process opens at lower density for APR HI. 

FigureEjshows the pressure-density (left) and the mass- 
central density (right) diagrams for all five EOSs. According 
to the left panel, all EOSs are similar. The P{p) curves in- 
tersect in one point, p ~ 8.5 x 10^* g cm~'^. However, small 
differences in the EOSs lead to noticeable differences in neu- 
tron star configurations, particularly, in the mass-density re- 
lations M{pc) (the right panel). Corresponding proton frac- 
tions and direct Urea thresholds are also different (Table 2 
and Fig. 1). 



4.2 Beyond the direct Urea threshold 

In contrast to the cooling scenarios in Papers I and II we 
use the EOS which opens the direct Urea process in massive 
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Table 2. Neutron star models for the EOSs DH, APR I, II, III, and APR 



Model 


Main 


DH 


APR 


APR I 


APR II 


APR III 




parameters 






7 = 0.6 


0.575 


0.643 


Maximum 




2.05 


2.20 


1.923 


1.919 


1.929 


mass 




2.86 


2.74 


2.759 


2.774 


2.731 


model 


R km 


9.99 


9.99 


10.32 


10.28 


10.39 


Direct Urea 


Md/M0 




2.01 


1.829 


1.891 


1.685 


threshold 






1.56 


1.680 


2.084 


1.275 


model 


R km 




10.95 


11.27 


10.83 


11.83 




2.2 



14.4 14.6 14.8 15 15.2 15.4 
Ig p (g cm-3) 




15 15.2 
Ig (g cm-3) 



15.4 



Figure 3. Pressu re — density (left) and mass - central density (right) diagrams for five EOSs. The APR EOS was suggested by 
Akmal et all il998h. The EOSs A PR I, II, and III refer to diff erent values of 7 (se e Table 2) in the parameterization proposed by 
Heiselberg fc Hiorth-JensenI Jl99Sft . The DH EOS was derived bv lDouchin fc Haense] 12001) . 



neutron sta rs. Moreover, contrary to our older cooling sce- 
narios (e.g.. lKaminker et alJl2002l[lYakovlev fc Pethickli2004l 
and references therein) the direct Urea process is now open 
at densities at which proton superfluidity dies out (Fig. 0, 
in agreement with some recent calculations of proton pair- 
ing in dense mat ter (e.g.. lTakatsuka fc Tamagakil997l . l2004l 
IZuo et alJliool . 

Figure0]shows a family of cooling curves similar to that 
in Fig. 121 It is calculated for the same models of proton and 
neutron pairing (p and nt, Fig. 1) but adopting the APR 
I EOS. The curves show the cooling of neutron stars with 
masses from 1.01 to Mmax = 1-92 Mq. The density of cooling 
curves reflects the probability to detect such stars provided 
the mass distribution of stars is flat in the indicated mass 
range (and observational selection effects are absent). Stars 
with M ^ Md = 1.83 M0 would have temperatures nearly 
in the same range as in Fig.|l(see Section 3). By opening the 
direct Urea we get a new but rare population of very cold 
neutron stars (Md 55 M 5C Mmax). Taking into account the 
difflculties to detect such stars, the chances to observe them 
are really low. Notice that if the EOS would allow higher 



Mmax we could obtaiu a new dense bunch of curves for very 
cold stars with M « Mmax- 

The same effects are demonstrated in Fig. |K| It shows 
Ta°° versus M (left) or pc (right) for stars of the Vela pulsar 
age, t = 1.66 x 10'* yrs. We present the curves for super- 
fluid stars with the EOSs APR I, II, and III. The horizontal 
dotted lines give the observational limits of for the Vela 
pulsar. 

One can see a very sharp fall of in the narrow mass 
range from M = Md to M = Mu + O.OIM© for the APR I 
EOS and especially for APR II. This fall is naturally associ- 
ated with the strong enhancement of the neutrino luminosity 
by the direct Urea process. 

The case of ARP III is different. As seen in Fig. Q the 
threshold of opening the direct Urea process is placed now in 
a wide region near the maximum of the critical temperature 
7cnt(p). When the internal stellar temperature falls notice- 
ably below Tent (p) in some region of p > pD , the suppression 
of the direct Urea process by neutron superfluidity becomes 
important and smears out the dependence of on M. The 
similar effect of partial suppressio n of the direct Urea proces s 
by superfluidity was discussed bv lPage fc Applegatj il992l) . 
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Figure 4. Same as in Fig. El but for the EOS APR I. Ninety two cooling curves refer to neutron stars of masses from 1.01 to 
Mniax = 1.92 Mq (with the step O.OIMq). The sharp boundary between densely and rarely covered domains corresponds to the direct 
Urea threshold, Md = l.SSM© (Table 2). 
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Figure 5. Surface temperatures of neutron stars of the age of Vela pulsar, t = 1.66 X 10^ yr, versus stellar mass (left) and central 
density (right) for three EOSs (Table 2): APR I (solid lines), APR II (short dashes), and APR III (thin and thick long dashes), assuming 
either combined proton and neutron (p and nt) pairing in stellar cores (thick lines) or only proton pairing alone {No n-SF; thin lines). 
The vertical dot-and-dashed lines, denoted by I, II, and III, show the threshold masses Mjj (left) or densities po (right) for the direct 
Urea process. Horizontal dotted lines are observational limits on the Vela's surface temperature. 
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Evidently, the level of this suppression is regulated by the 
strength and density range of neutron pairing and by the 
threshold density of the direct Urea process. Were neutron 
superfluidity absent (the thin long-dashed line), the depen- 
dence would be as sharp as for the EOSs ARP I and II. 

The dependence of Tg°° on pc (the right panel of Fig. |3 
is similar to the dependence on M. Note only much smoother 
dependence of on pc in the lower parts of the curves; 
small variations of M < Mmax correspond to larger varia- 
tions of Pc (Fig. EJ- 

Figure |S| is similar to Figs. |5| and 2] and shows cooling 
curves calculated for the APR III EOS and the same models 
for proton and neutron pairing (p and nt). Neutron star 
masses range from 1.01 to Mmax = 1.93 Mq. A sharp lower 
boundary of the densely covered domain corresponds to the 
direct Urea threshold Mb = I.68M0. In contrast to Fig. 01 
we now obtain a representative (densely populated) family 
of middle-aged neutron stars with M > Mo and lgT!°° ~ 
5.6 — 5.7. They are much colder than stars with M < Mo but 
much hotter than those with M ~ A/max. They appear owing 
to the flattening of the dependence on M by neutron 
superfluidity (Fig. |3 . Were this cooling scenario realized in 
the nature, the chances to detect such stars would be higher 
than for the scenario in Fig. 2] 

The comparison of Figs. |1] and El shows a slight increase 
of the gap between the lower boundary T^{t) of densely 
populated family of neutron stars with M < Mo and the 
observed surface temperature of the Vela pulsar. This dif- 
ference is also illustrated in Fig. |^ the curves (M) and 
Ta°°(pc) merge in nearly the same slowly decreasing function 
before the direct Urea thresholds (M < Mo, pc < po)- 



5 DISCUSSION AND CONCLUSIONS 

We have extended the scenario of neutron star cooling pro- 
posed in Papers I and II to the class of EOSs which allow the 
operation of the powerful direct Urea process in the cores of 
massive stars. It widens possible interpretations of observa- 
tional data on thermal emission from isolated middle-aged 
neutron stars using simplest models of neutron star cores 
composed of nucleons, electrons, and muons. 

The main common feature of our present scenario and 
the scenarios of Papers I and II is that the enhanced neu- 
trino emission, required for the interpretation of observa- 
tions of coldest neutron stars, is produced by Cooper pairing 
of neutrons. This puts stringent constraints (Paper I) on the 
density dependence of the critical temperature Tcnt(p) for 
triplet-state neutron pairing. Tuning our phenomenological 
model of rcnt(p) we have obtained a noticeable dependence 
of cooling on neutron star mass M. 

In the scenario of Paper I we have also needed strong 
proton pairing in the outer cores of neutron stars to explain 
the observations of neutron stars hottest for their age. In Pa- 
per II we have extended that scenario taking into account 
a possible presence of accreted envelopes and singlet-state 
pairing of neutrons in the stellar crust. We have shown that 
the presence of accreted envelopes allows us to take weaker 
proton pairing. Note, that models of weaker pairing are con- 
sistent with recent microscopic ca lculatio ns of proton criti- 
cal t emperatures ijZuo et alJ I20fl 4. Takat ~ka fc Tamaeakil 
I2OO4I) . although some other calculations predict much 



stronger proton pai ring ('e.g.. [Lo rnbardo fc Schulz3 1200 ll : 
also see references in lYakovlev et al.lll999l . and a recent pa- 
per bv lTanigawa et al.ll20o3i ! 

In this paper, for simplicity, we use the model of 
strong proton pairing in the stellar core. However, follow- 
ing Paper II, we could assume the presence of accreted 
envelopes, which would again allow us to weaken proton 
pairing. Adopting moderately strong proton pairing (with 
r™''" > 10® K) we could interpret the observations of neu- 
tron stars warmest for their age. 

The main feature of this paper is the consideration of 
EOSs which open the direct Urea process in central layers 
(P <^ 10^^ g cm~^) of most massive star s. We have employed 
the paramet erization of the EOS o f lAkmal et alj l| l998h 
proposed by iHeiselberg fc Hiorth-JensenI Jl999l l200(i) . In 
combination with properly chosen phenomenological mod- 
els for strong proton pairing {T^'^^ ^ 5 x 10® K at 
p < 8 X 10 g cm" ) and moderate neutron pairing (r^.^^ ~ 
6 X 10* K at p > 10^^ g cm~^) it allows us to predict new 
types of cooling neutron stars, in addition to those given by 
the cooling scenario without any direct Urea process (Pa- 
pers I and II, Sect. 3). They are rapidly cooling neutron stars 
with masses higher than the direct Urea threshold mass Mo ■ 
We have shown (Section 4) that the family of these stars can 
be either rarely populated (Fig. jlj and thus difficult to dis- 
cover, or rather distinct and populated (Fig. |SJ , depending 
on the choice of a model EOS of dense matter and a model 
for neutron superfluidity. Neutron stars from the distinct 
class could be cold but much warmer than massive stars 
with the fully open direct Urea process. 

There is no clear evidence for the existence of such cold 
stars at the moment but it may appear in the future. For in- 
stance, low upper limits on thermal emission fr om hypotheti- 
cal neutron stars in a few supernova remnants (iKapl an et alJ 
[2004} can be treated as an indirect indication of the exis- 
tence of such cold stars. In this context the new scenario 
would allow one to explain the existence of these cold ob- 
jects (if discovered) in the frame of standard and simple 
physics of dense matter (without assuming the presence of 
hyperons and/or exotic forms of matter, such as pion/kaon 
condensates or quarks). 

As has already been emphasized, we need neutron pair- 
ing to explain the current observa tions of stars co l dest fo r 
their age. However, according to 'Ousak ov et alJ (l2004bl) . 
cooling curves are not too sensitive to exchanging T'cp(p) ^ 
rcnt(p) in neutron star cores. Therefore, we would also be 
able to explain the data assuming strong neutron pairing 
and moderate proton pairing in stellar cores. 

Finally, let us note that the same physics of neu- 
tron star interiors, which is tested by observations of iso- 
lated (cooling) neutron stars, can also be tested by ob- 
servations of accreting neutron stars in soft X-ray tran- 
sients (e.g., Yakovlev, Levenfish & Haensel 2003) basing 
on the hypothesis of deep crus tal heating of such stars 
iBrown. Bildsten. Rutledgell988^ bv pycnonu clear reactions 
in accreted matter faaensel fc Zduniklll99ot) . The observa- 
tions of soft X-ray transients in quiescent states indicate 
(Yakovlev, Levenfish & Gnedin 2005) the existence of rather 
cold neutron stars (first of all, SAX J1808. 4-3658) inconsis- 
tent with the minimal model of cooling neutron stars. Al- 
though these observa tions are currently inconclusive (e.g., 
lYakovlev et"ai]l2005l) . if confirmed in future observations. 
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Ig t (yrs) 

Figure 6. Same as in Figs. 1^ and l^but for the APR III EOS. Ninety three thin solid curves refer to neutron stars of masses from 1.01 
to Afniax = 1.93 M0 (with the step 0.01 Mq). The lower boundary of the densely covered domain refers to the direct Urea threshold, 
Md = 1.68Mq. 



they could give an evidence in favor of the cooling scenario 
proposed in this paper. 
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